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Q ; ABSTRACT 

From a search of a portion of the sky covered by the SDSS and UKIDSS 
databases, we have located 2 galaxies at z ~ 0.5 that have properties similar 
to those of the luminous passive compact galaxies found at z ~ 2.5. From 
Keck moderate-resolution spectroscopy and laser-guided adaptive-optics imaging 
of these galaxies, we can begin to put together a more detailed picture of what 
their high-redshift counterparts might be like. Spectral-synthesis models that 
fit the u to K photometry also seem to give good fits to the spectral features, 
i— 1 1 From these models, we estimate masses in the range of 3-4 x 1O 11 M for both 

\ galaxies. Under the assumption that these are spheroidal galaxies, our velocity 

dispersions give estimated masses about a factor of 3 smaller. However, our high- 
t^j- \ resolution imaging data indicate that these galaxies are not normal spheroids, 

and the interpretation of the kinematic data depends critically on the actual 
\ morphologies and the nature of the stellar orbits. While recent suggestions that 

i \ the population of high-redshift compact galaxies is present locally as the inner 

q \ regions of local massive elliptical galaxies are quite plausible, the peak mass 

surface densities of the two galaxies we discuss here appear to be up to a factor 
of 10 higher than those of the highest density local ellipticals, assuming that 
our photometric masses are roughly correct. It thus seems possible that some 
dynamical "puffing-up" of the high-redshift galaxies might still be required in 
this scenario. 

Subject headings: galaxies: high-redshift — galaxies: evolution — galaxies: structure- 
galaxies: stellar content 



1 Somc of the data presented herein were obtained at the W.M. Keck Observatory, which is operated 
as a scientific partnership among the California Institute of Technology, the University of California and 
the National Aeronautics and Space Administration. The Observatory was made possible by the generous 
financial support of the W.M. Keck Foundation. 
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Introduction 



Luminous passive galaxies at z > 2 general ly are very compact, with R e < 2 kpc (e.g. 



van Dokkum et al.l 120081 ; iDamjanov et al.l l2009h. Many of these have nearly expon ential 



profiles and appear to be disk-like (jStockton et al.l 12004 ; IStockton fc McGrathl 120071 ) , but 
at least some of the most compact ones, typi cally with R P < 1 kpc, s eem to have surface- 
brightness profiles closer to an r 1//4 law (e.g., Ivan Dokkum et al.l 120081 ) . The most extreme 
example that we know of is a galaxy with an apparent mass of ~ 3 x 10 11 M Q and an 
R e = 440 pc at z = 2.48 (Stockton et al., in preparation). Galaxies such as these are absent, 
or at least extremely rare, in the present-day universe, and much of the recent discussion has 
centered on either the reliability of the size estimates at high redshifts or the evolutionary 
path such galaxies might h ave taken and how their stars have been incorporated into galaxies 



in th e local universe (e.g. , iDavel 120081 ; ICimatti et al.l 120081 ; Fan et al.l 120081 ; iBezanson et al. 



20091 ; iHopkins et al J 12009( 1 . 



It is extremely difficult to study these luminous, compact galaxies at z ~ 2.5 beyond 
obtaining their photometric SEDs, which confirms their approximate redshifts, and high- 
resolution images (HST or adaptive optics), which can indicate their gross morphologies. 
Photometric masses can be estimated from model fits to the SEDs, but these are uncertain 
because of partial degeneraci es among age, metallicity, and r eddening, as well as basic uncer- 
tainties in the models (e.g., iMuzzin et al.l 120091 ) . Recently, Ivan Dokkum et al.l (120091 ) have 
claimed a very high velocity dispersion of 510 km s _1 for a compact galaxy at z — 2.2 from 
29 hours of observation with GNIRS on Gemini South, but this observation only empha- 
sizes the heroic efforts that are needed with current facilities to obtain such results, which, 
even so, are not immune to worries that unrecognized systematic effects may jeopardize any 
conclusions reached. It would be extremely useful to find similar galaxies at lower redshifts, 
where they could be studied in greater detail. 



There have been some recent efforts to identify such galaxies. iTaylor et al.l (120091 ) 
searched the Sloan Digital Sky Survey (SDSS) database for early-type compact galaxies 
in the redshift range 0.066 < z < 0.12, using both SDSS spectra and photometric redshifts. 
They found a number of galaxies with 1.1 < R e < 1.5, but these all had indicated masses 
of ~ 5 x 10 10 Mq, so they are n ot really comparable to the massive compact galaxies found 
at high redshifts. ITaylor et al.l conclude that such galaxies must be extremely rare at the 
present epoch and that their size evoluti on cannot be a resu l t of a stochastic mechanism such 
as major mergers. On the other hand, IValentinuzzi et al.l (120091 ) claim to find substantial 
numbers of massive, old, compact galaxi es in nearby X-ray-s elected clusters, including a few 
with M > 10 11 M Q and R e < 1.5 kpc. IValentinuzzi et al.l use a circularized value for R e , 
rather than the more common major axis value (which we use throughout this Letter), but 
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it is not clear how much this affects their conclusions. 

We present here some initial results from a search at z ~ 0.5 of a portion of the area 
common to the SDSS and the UKIRT Infrared Sky Survey (UKIDSS) surveys, describing in 
some detail results we have obtained for two of the galaxies that have properties within the 
range of the passive compact galaxies found at high redshifts. Throughout, we assume a flat 
cosmology with H = 71 km s _1 Mpc -1 and Q m = 0.27. 



2. Object Selection and Observations 

We selected objects from the area common to the SDSS DR7 and the UKIDSS DR5plus 
for RAs between 21 h and 2 h , a total area of ~ 200 deg 2 . In this region, we searched the 
combined UKIDSS/SDSS database for objects (1) that had SDSS mj (model) — m^PSF < 0.3 
(this criterion ensures that the best-fitting model for the object is close to a PSF profile), 
(2) that had colors expected for old stellar populations at a given redshift, and (3) that 
had magnitudes that would place them at or above ~ 2L* at that redshift. We made 
separate searches at redshifts ranging from 0.40 to 0.60, stepping in intervals of 0.05 in z. 
To pass initial selection, an object had to match an old-stellar population SED in 7 colors 
constructed from adjacent filters in the SDSS and UKIDSS surveys (no colors crossing the 
surveys were used), where the fitting range was ±3 times the typical uncertainty, estimated 
separately for each color. The only serious contaminating objects we have identified are 
certain varieties of carbon stars, which tend to mimic old-stellar-population colors at z ~ 0.4 
surprisingly well over the 9 bands we have from the UKIDSS/SDSS databases. We made a 
final prioritization of objects for furt her observation by evaluating the galaxies' profiles in the 



SDSS i-band images with galfit (IPeng et al.l 120021 ) and eliminating those with apparent 



R e > 1.5 pixels. The final sample comprised 14 objects. The two galaxies we discuss 
here, SDSS J232949.60+151106.3 (hereinafter SDSSJ2329) and SDSS J011436.33-011438.1 
(hereinafter SDSSJ0114), are both classified as "stars" in the SDSS database. 

We carried out spectroscopy of 11 of the 14 candidate galaxies on 2 009 Aug 22 UT and 



2009 Sep 16 UT with the Low- Resolution Imaging Spectrograph (LRIS; lOke et al.lll995l ) on 
the Keck I telescope. The spectra were obtained with the 600/5000 grating and a l"-wide 
slit. The detector consisted of 2 LBL 2048 x 4096 CCDs, binned 2x2. This configuration 
gives a resolution of ~ 5 A. Most of the objects were, in fact, galaxies at redshifts ranging 
from 0.4 to 0.8. Three were carbon stars, and 1 was a peculiar QSO. 



We also obtained if-band images of 6 of these galaxies on the night of 2009 Sep 15 UT 
with the Keck II laser-guide-star adaptive-optics system (LGSAO; IWizinowich et al.l 120061 ) 
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and the NIRC2 camera. Total exposures for each of the two objects discussed here were 18 
min at the 40 mas pixel scale. Before registering and co-adding the individual exposures, 
we corrected them for distortion onto a 2x oversampled grid, using the drizzle algorithm 
fiFruchter k Hookll2002h with a pixel reduction scale of 0.8. The PSF FWHM was ~ 80 mas, 
limited by the 40 mas scale. Both of our fields have bright (r < 18.1) stars within ~ 16" 
of the galaxies that can be used for PSF determination, and our exposures were kept short 
enough to avoid saturating on these stars. Observations of a globular cluster field on the 
same night indicate that, for our purposes, PSF variations over this distance were negligible. 
Comparison of previous similar AO observations onaz~ 2.5 galaxy with HST NICMOS 
obs ervations of the same galaxy gave essentially identical values for the Sersic n index and 
R e (jStockton et al.ll2008l ). so we can have reasonable confidence in the detailed results from 
our AO imaging. 



Results 



3.1. Spectral-Energy Distributions and Stellar Populations 



We first show in Fig. CD the spectral-energy distributions (SEDs) for SDSSJ2329 and 
SDSSJ0114, based on their SDSS/UKIDSS magnitudes. Because o f the possibility tha t the 
zeropoints of the two surveys may differ at the few percent level (IHewett et al.l 120061 ). we 
have imposed a minimum uncertainty of 5% on the magnitudes. We have used a modified 
version of the photometric redshift code Hyper-z (IBolzonella et al.l 120001 ) to fit SEDs to 
the magnitudes, where we have constrained the redshifts to be those we have determined 
from the spectroscopy described below. We explored a range of models, including both 
instantaneous-bur st models and model s with exponentially decreasing star- formation rates, 
all with a range of ICalzetti et al.l (120001 ) reddening and metallicities of 0.4, 1.0, and 2.5 solar. 
For both galaxies, the best fitting SED was from an instantaneous-burst model with no 
reddening. For SDSSJ2329, the model had an age of 5.0 Gyr and a metallicity of 0.4 solar, 
while for SDSSJ0114, the model had an age of 4.75 Gyr and solar metallicity. The indicated 
photometric masses, determined from the current mass in stars from the HYPER-Z model 
fits, are 3.9 x 1O U M for SDSSJ2329, and 3.4 x 10 n M Q for SDSSJ0114. 

Figure [2] shows the LRIS spectra of SDSSJ2329 and SDSSJ0114, compared with the 
models that give the best fits to the photometry, as shown in Fig. [TJ Here we have simply 
scaled the models to the spectra (to account for slit losses). The redshifts for these two galax- 
ies are fortuitously virtually the same, with SDSSJ2329 having z = 0.490 and SDSSJ0114 
having z = 0.489. 
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Fig. 1. — Hyper-z best-fit SEDs, constrained by the spectroscopic redshifts for the galaxies 
and selecting among Charlot-Bruzual (2007; privat e commu n icatio n) instantaneous burst 
models with 0.4, 1.0, and 2.5 solar metallicities and IChabrierl (120031 ) initial mass functions. 
SDSSJ2329 and SDSSJ0114 were fit best by 0.4 solar and 1.0 solar models, respectively. 
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Fig. 2. — Spectra of SDSSJ2329 and SDSSJ0114 (black traces) superposed on the spectral 
synthesis models that best fit the SDSS/UKIDSS photometry (red [gray] traces), as shown 
in Fig. [TJ 

3.2. Morphologies 

From our LGSAO imaging, we were able to explore the morphologies and radial- surf ace- 
brig htness profiles of SPSS J2329 and SDSSJ0114. Using the model-fitting application gal- 



fit ( iPeng et al.ll2002l ) and PSFs from nearby stars in the field, we first fit single Sersic models 



to the galaxy images, as shown in Fig. [3] (second panels). These fits give Sersic R e = 0.92 
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Fig. 3.— Keck LGSAO images and model fits for SDSSJ2329 and SDSSJ0114. For each 
galaxy, the panels from left to right show the original LGSAO image, the residual from the 
best-fit single-Sersic model, the residual from the best-fit double Sersic model, and the double 
Sersic model itself, at the same contrast as the original image. The insets in the model panels 
show the models without convolution with the PSFs, at slightly lower contrast. The Sersic 
index n and the major axis R e are given below the corresponding model panels. 

kpc for SDSSJ2329 and R e = 1.37 kpc for SDSSJ0114, in the range of the compact passive 
galaxies found at high redshifts, although not as small as the most extreme examples, which 
have R e ~ 0.5 kpc. These single-component fits also leave systematic residuals with peak 
values ~ 2% of the peaks of the original images. We therefore explored 2-component fits, as 
shown in Fig. [3] (third and fourth panels). These, as expected, show much smaller residuals; 
the main point, however, is that the residuals show little or no systematic structure. 

An interesting feature of these fits is that single Sersic models give Sersic indices n > 4, 
typical of spheroidal galaxies. At high redshifts, the small but systematic residuals left by 
these fits would fade into the noise. For the two-component Sersic fits, which give much 
more satisfactory residuals, the Sersic indices for both components are closer to exponential 
profiles. Furthermore, for the more extended components of our 2-component Sersic fits, the 
projected axial ratios are fairly small: b/a = 0.34 for SDSSJ2329 and 0.40 for SDSSJ0114. We 
have tried replacing the extremely compact central components with both a PSF component 
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Fig. A.— Radial-surface-brightness profiles for SDSSJ2329 and SDSSJ0114, from the LGSAO 
images shown in Fig. 3. The solid red (gray) lines show the best Sersic 2-component fit, with 
the green and blue (light gray) curves showing the individual components. The lower panels 
show the differences between the annulus photometry and the model fits. 



and a small r 1//4 -law bulge, but the fits were significantly worse for both of these models. 

Figure H] shows the radial-surface-brightness profiles along the major axes, derived from 
elliptical annuli sampling of the images, along with the best 2-compone nt model fits. F igure^] 
shows the corresponding mass surface- density profiles (assuming a IChabrierl 120031 IMF), 
together with the "upper-envelope" profile of iHopkins et al.l (120091 ). which represents an 
approximate upper limit to the mass surface density at each radius for any galaxy from an 
ensemble of massive elliptical galaxies in the local universe. 



3.3. Velocity Dispersions 



We have estimated the vel ocity dispersions fro m our spectra via the direct method, using 
spectra of G giant stars from IValdes et al.l (120041 ) as templates. After re-binning all of the 
spectra to velocity units, we broadened the template spectra with an appropriate Gaussian 
to bring them to the same resolution as our de-redshifted spectra. For each galaxy, a grid 
of template spectra was prepared with a range of Gaussian broadening parameters, and the 
scaled template spectra were subtracted from the galaxy spectra. After subtraction of a 
low-order spline to take care of any local mismatches between the template and the galaxy 
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Fig. 5. — Mass-surface-density profiles for SDSSJ2329 and SDSSJ0114, assuming that the 
photometric masses derived from the spectral-synthesis model fits are correct. The radial 
coordinate R is along the semi-major axis, and the surface densities are measured in elliptical 
annuli from the 2-component Sersic models, without convolution with the PSFs, which should 
closely approximate the true global morphologies of the galaxies beyond a radius of ~ 300 pc 
(c orresponding to R 1 ^ ~ 0.75 kpc 1 / 4 ). The dotted line shows the "upper-envelope" profile 



of iHopkins et al.l (120091 ) . which gives the approximate upper limit to the mass density at a 
given physical radius from a sample of massive elliptical galaxies at low redshifts. The range 
shown is 0.11 < R < 8.8 kpc. 

spectra, the \ 2 statistic was calculated for each broadening parameter. The minimum x 2 
value was found by fitting a curve to \ 2 values in the vicinity of the minimum. In order 
to estimate the errors, we did Monte-Carlo simulations, creating 100 mock spectra that 
simulated both the random and correlated noise properties of the actual spectra for each 
galaxy. 

This procedure gives a = 291 ± 16 km s _1 for SDSSJ2329 and a = 254 ±21 km s _1 
for SDSSJ0114. For spheroids, velocity dispersions can be translated to estimates of the 
dynamical mass by Md yn = kR e a 2 /G, where k is a parameter that depends on the properties 
of the stellar velocity distribution function for the galaxy, as well as other fact ors, such as the 



regio n of the galaxy covered by the spectrum. It is often taken to be ~ 5 (e.g-. lBernardi et al. 



20091 1. With this assumption, we find that M dyn = 0.91 x 1O U M for SDSSJ2329 and 1.0 x 



10 M & for SDSSJ0114, where we have used the single-Sersic estimates for R e . These values 
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are about 30% of the masses estimated from the stellar population fits. The discrepancy 
could be attributed to problems with the assumption of a standard IMF or concerns about 
the value used for R e (although, since the latter enters linearly, it is unlikely to make up the 
entire difference). 

However, it is also possible that it is simply wrong to assume that a procedure normally 
applied to ellipticals is wholly relevant in these cases. Both the Sersic indices of the best- 
fitting models and the axial ratios of more extended components fall outside the range of 
normal spheroids, so our "velocity dispersion" could actually be dominated by bulk velocities 
of one kind or another. We have modeled rotation curves, under the assumption that both 
components of each of these galaxies are rotating disks, that the mass-to-light ratio is con- 
stant, and that the disks can be approximated by flattened spheroids. This last assumption 
simplifies the calculation and is a close enough approximation to the exp onential disk case to 



show the general form of the expected profile fairly accurately (see, e.g.. lBinney fc Tremaine 



19871 ). These models show that even for total masses < 0.1 of the photometric mass we have 
found, we would expect to see double lines at the resolution of our integrated spectra, which 
are clearly not present. Another possibility is that the galaxies are prolate, with mostly 
radial orbits aligned with the long axis. Further work, including resolved spectroscopy, is 
needed to explore this issue. 



4. Discussion 



The two galaxies we discuss here show many of the characteristics of the luminous 
compact passive galaxies found at high redshift, and their relative proximity allows us to 
explore them in much greater detail than is possible for the high-redshift cases. However, they 
cannot actually be survivors of the population found at z > 2, since their stellar population 
ages indicate that most of the stars in them were formed at z ~ 1.8. Nevertheless, it is still 
quite possible that they have had similar formation histories to those of the high-redshift 
population, although perhaps in less dense environments, where formation processes might 
have been delayed. Assuming that we can take SDSSJ2329 and SDSSJ0114 as mirroring the 
properties of the compact high-redshift galaxies, what implications do they suggest? 



Our surface-brightness profiles shown in Fig. H] cover a range of over 6 magnitudes over 
a semi-major axis of ~ 6 kpc. The fact that we can probe these galaxies to much lower 
surface-brightness levels than we can their high-redshift counterparts shows that, in at 
least these cases, their compactness is unlikely to be explaine d by suggestions tha t we 



are simply missing low-surface-brightness outer regions (e.g.. iHopkins et al.ll2009l ). 
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• The necessity of two-component fits to avoid systematic residuals indicates that these 
galaxies have already undergone considerable dynamical evolution. They are not as 
simple systems as might have been assumed from the available data on their high- 
redshift counterparts. 

• The best fitting models are superpositions of components with surface-brightness pro- 
files close to exponentials, even though the single-component models have profiles close 
to r 1 / 4 laws. In addition, the more extended components of the reconstructed models 
have high axial ratios, as seen in the insets to Fig. [3j These conclusions line up with 
evidence that luminous passive galaxies at high redshifts often have exponential pro- 
files and an apparent disk- like appearance (IStockton et al.ll2004l ; IStockton fc McGrath 



20071 ; IStockton et al.ll2008l ). However, the lack kinematic indications of a disk in our 



integrated spectra suggest that these galaxies may not, in fact, be rotationally sup- 
ported. Prolate morphologies with strongly anisotropic stellar velocity fields remain a 
possibility. 

The two-component models that best fit the luminosity profiles indicate that slightly 
more than half of the light from SDSSJ2329 and SDSSJ0114 comes from structures 
that have R e ~ 2 kpc, still small, but considerably larger than the R e found from 
single-component models. On the other hand, > 1/3 of the ligh t comes from ext remely 
compact structures with R e < 250 pc. Assuming a standard IChabrierl (120031 ) initial 
mass function (IMF), our spectral synthesis model fits in Fig. [1] imply masses of > 
1O 11 M for these extremely compact cores. 



Bezanson et al.l (120091 ) and iHopkins et al.l (120091 ) present plausible and convincing ar- 



guments that the luminous passive compact galaxies at high redshifts are locally present as 
the core regions of the most massive present-day elliptical galaxies, which have presumably 
acquired their exten ded envelopes through mostly dry mergers with more recently formed 
and fluffier galaxies. IHopkins et al.l (120091 ) define an "upper envelope" to the surface mass 
density profiles of massive local ellipticals and show that the apparent mass surface density 
of the high-redshift compact galaxies rarely rises above this profile at the same physical radii. 
We have plotted this upper-envelope curve in Fig. [5j for comparison with the mass-surface- 
density curves for our best-fit photometric mass models for SDSS J2329 and SDSSJ0114. The 
surface densities for both models are a factor of ~ 3 higher than the IHopkins et al.l (120091 ) 
upper envelope even out to R = 2 kpc. This factor rises to ~ 10 in the inner regions. If 
this difference can be trusted (and if it is representative of the high-redshift counterparts to 
these galaxies), it indicates that there still needs to be some amount of "puffing up" of these 
objects to bring their central surface densities into line with those of massive ellipticals at 
low redshift. Of course, our mass-surface-density profiles might be overestimated if our total 
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mass estimates are too large because (1) we have grossly overestimated the ages of the stellar 
populations (we believe t his to be u nlikel y) or (2) the IMFs for the stellar populations are 
more top-heavy than the IChabrierl (120031 ) IMFs we have assumed. If the actual masses of 
these galaxies were determined to be more in line with our velocity dispersion measurements, 
this latter possibility would provide the most likely explanation. Future observations of these 
and other similar galaxies at moderate redshifts should give us insight into these issues. 



We thank the anonymous referee for a careful and thoughtful reading of the original 
version of this Letter and for offering numerous suggestions to improve both its substance 
and presentation. We also th ank Phil Hopkin s for h is comments on the earlier version. The 
UKIDSS project is defined in lLawrence et all toOlh . UKIDSS uses the UKIR T Wide Field 



Came ra (WFCAM; ICasali et al.l 120071 ) and a photometric system des cribed inlHewett et al. 



( 2006 ). The pipeline processing and science archive are described in Hambly et al. (j2008l ). 
Funding for the SDSS has been provided by the Alfred P. Sloan Foundation, the Participating 
Institutions, the National Science Foundation, the U.S. Department of Energy, the National 
Aeronautics and Space Administration, the Japanese Monbukagakusho, the Max Planck 
Society, and the Higher Education Funding Council for England. The SDSS Web Site is 
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